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ABSTRACT 

A recently developed micromechanical theory for the thermoelastic response of function- 
ally graded composites with nonuniform fiber spacing in the through-thickness direction is 
further extended to enable analysis of material architectures characterized by arbitrarily nonuni- 
form fiber spacing in two directions. In contrast to currently employed micromechanical 
approaches applied to functionally graded materials, which decouple the local and global effects 
by assuming the existence of a representative volume element at every point within the compo- 
site, the new theory explicitly couples the local and global effects. The analytical development is 
based on volumetric averaging of the various field quantities, together with imposition of boun- 
dary and interfacial conditions in an average sense. Results are presented that illustrate the capa- 
bility of the derived theory to capture local stress gradients at the free edge of a laminated com- 
posite plate due to the application of a uniform temperature change. It is further shown that it is 
posssible to reduce the magnitude of these stress concentrations by a proper management of the 
microstructure of the composite plies near the free edge. Thus by an appropriate tailoring of the 
microstructure it is possible to reduce or prevent the likelihood of delamination at free edges of 
standard composite laminates. 


NOMENCLATURE 

indices used to identify the cell ip,q,r) 
number of cells in the x 2 and x 3 directions, respectively 
indices used to identify the subcell (afty) 
dimensions of the subcell (aPy) in the (p,q,r)lh unit cell 
volume of the subcell (afty) in the (p,q,r)\h unit cell 

’On leave from Tel- Aviv University, Ramat-Aviv 69978, Israel. 
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local subcell coordinates 

coefficients of heat conductivity of the material in the subcell (apy) 
temperature field in the subcell (aPy) 

temperature at the center of the subcell (aPy) when / = m = n = 0; 
coefficients associated with higher-order terms in the temperature field 
expansion within the subcell (aPy) for other values of l,m,n 

components of the heat flux vector in the subcell (aPy) 

average values of the subcell heat flux component when 

l=m=n= 0; higher-order heat fluxes for other values of l,m,n 

surface integrals of subcell interfacial heat fluxes 

displacement components in the subcell (aPy) 

x i displacement components at the center of the subcell (apy) when 
l = m, = n = 0; coefficients associated with higher-order terms in the 
displacement field expansion within the subcell (aPy) for other values of 
l,m,n. 

local strain components in the subcell (apy) 

local stress components in the subcell (apy) 

elements of the stiffness tensor of the material in the subcell (aPy) 

elements of the thermal tensor of the material in the subcell (aPy) 

average values of the subcell stress components when 

/ = m = n = 0; higher-order stress components for other values of l,m,n 

surface integrals of the subcell interfacial stresses at x j =±d a / 2 

surface integrals of the subcell interfacial stresses at xi = ±fcjf V2 

surface integrals of the subcell interfacial stresses at x-i, = ±/y V2 
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1.0 INTRODUCTION 

Functionally graded materials (FGMs) are a new generation of composite materials in 
which the microstructural details are spatially varied through nonuniform distribution of the 
reinforcement phase, by using reinforcement with different properties, sizes and shapes, as well 
as by interchanging the roles of reinforcement and matrix phases in a continuous manner. The 
result is a microstructure that produces continuously changing thermal and mechanical proper- 
ties at the microscopic or continuum level. 

The use of functionally graded materials in applications involving severe thermal gradients 
is quickly gaining acceptance in the composite mechanics community and the aerospace and air- 
craft industry. This is particularly true in Japan and Europe, where the concept of FGMs was 
conceived. The current approach employed by the Japanese and European researchers in analyz- 
ing the response of FGMs to thermal gradients is the standard micromechanics approach based 
on the concept of a representative volume element (RVE) assumed to be definable at each point 
within the heterogeneous material (cf., Wakashima and Tsukamoto, 1990; Fukushima, 1992). 
This assumption, however, neglects the possibility of coupling between local and global effects, 
thus leading to potentially erroneous results in the presence of macroscopically nonuniform 
material properties and large field variable gradients. This is particularly true when the tempera- 
ture gradient is large with respect to the dimension of the inclusion phase, the characteristic 
dimension of the inclusion phase is large relative to the global dimensions of the composite, and 
the number of uniformly or nonuniformly distributed inclusions is relatively small (Aboudi et 
al., 1993). Perhaps the most important objection to using the standard RVE-based micromechan- 
ics approach in the analysis of FGMs is the lack of a theoretical basis for the definition of an 
RVE, which clearly cannot be unique in the presence of continuously changing properties 
due to nonuniform inclusion spacing. 

As a result of the limitations and shortcomings of the standard micromechanics approach, a 
new higher order micromechanical theory for functionally graded materials, HOTFGM, that 
explicitly couples the local and global effects, has been developed (Aboudi et al., 1993; Aboudi 
et al., 1994a,b). The results obtained thus far have demonstrated that the theory is an accurate, 
efficient and viable tool in the analysis of functionally graded materials and design of function- 
ally graded architectures in metal matrix composites. These results include verification of the 
accuracy of HOTFGM using the finite-element method (Pindera and Dunn, 1994), and the 
assessment of the applicability of the uncoupled micromechanics approach for the analysis of 
functionally graded materials (Pindera et al., 1994a, b). In particular, comparison of results 
obtained using the standard micromechanics approach with those of HOTFGM has demon- 
strated the need for a theory which explicitly takes into account the micro-macrostructural 
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coupling effects, thus justifying the development of the coupled higher-order theory. 

HOTFGM is a recently constructed theory that continues to evolve. The original formula- 
tion has been developed in the Cartesian coordinate system, and was intended for the analysis of 
functionally graded plates subjected to a temperature gradient across the plate’s thickness that 
coincides with the direction along which the microstructure is graded. The most recent develop- 
ments of the Cartesian-based theory include incorporation of two inelastic constitutive models 
for the response of metallic matrices (Aboudi et al., 1994c) and extension of the theoretical 
framework to include generalized plane strain loading situations in order to facilitate modeling 
of actual functionally graded structural components (Aboudi, et al., 1995). 

In this paper we present a further extension of HOTFGM that involves development of a 
two-dimensional framework to enable modeling of composites functionally graded in two direc- 
tions. The analytical approach in the two-dimensional theory, as in the one-dimensional version, 
is based on volumetric averaging of the various field quantities together with the imposition of 
boundary and interfacial continuity conditions in an average sense. The previous restriction of 
periodicity in two orthogonal directions however, is presently abondoned thus allowing arbitrary 
distribution of one or more reinforcement phases in one plane. This leads to a significant gen- 
eralization of the theory. As a result, composites with finite dimensions along the functionally 
graded directions can be analyzed. Figure 1 illustrates the types of internal architectures that can 
be analyzed with this new two-dimensional version of HOTFGM. These architectures include 
rows of aligned inclusions (or continuous fibers) with variable spacing in the functionally graded 
x 2 and *3 directions and regular spacing in the periodic x x direction. Figure la. Alternatively, 
completely random inclusion (or fiber) architectures in the x 2 - x 3 plane can also be admitted. 
Figure lb. At present, the two-dimensional version of the theory, herein called HOTFGM-2D, 
is limited to the analysis of functionally graded composites in the linearly elastic range. 

This theory is subsequently employed to study the free-edge problem in a symmetrically 
laminated B/Ep-Ti composite plate subjected to a uniform temperature change. The capability of 
the theory to capture large stress gradients near a geometric discontinuity such as the free edge is 
established upon comparison with finite-element analysis carried out by Herakovich (1976) 
using homogenized properties for the B/Ep plies. Subsequent incorporation of the actual micros- 
tructure of the B/Ep plies in the HOTFGM-2D analysis of the free-edge stress fields demon- 
strates the limitations of the homogenized continuum approach in the presence of course micros- 
tructure and large stress gradients. Finally, the potential of using functionally graded fiber archi- 
tectures in reducing edge effects in laminated MMC plates is demonstrated by investigating the 
effect of nonuniform fiber distributions in the B/Ep plies near the free edge. It should be noted 
that even though the utility of the theory is demonstrated herein for the special case of a 
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symmetric laminate under uniform temperature change, the theory naturally can be employed in 
more complicated situations with non-zero temperature gradients. 


2.0 ANALYTICAL MODEL 

HOTFGM-2D is based on the geometric model of a heterogeneous composite with a finite 
thickness H, and finite length L, that is infinite in the x x direction (see Figure 1). The loading 
applied to the boundaries of the composite in the x 2 - x 3 plane may involve an arbitrary tempera- 
ture distribution and mechanical effects represented by a combination of surface displacements 
and/or tractions. The composite is reinforced in the x 2 -x 3 plane by an arbitrary distribution of 
infinitely long fibers oriented along the x\ axis, or finite-length inclusions that are arranged in a 
periodic manner in direction of the x { axis. The heterogeneous composite is constructed using a 
basic building block (p,q,r ), Figure 2a, consisting of eight subcells designated by the triplet 
(aPy), Figure 2b. Each index a, p, y takes on the values 1 or 2 which indicate the relative position 
of the given subcell along the jcj, x 2 and x 3 axis, respectively. The dimensions of the unit cell 
along the xi axis, d u d 2 , are fixed for the given configuration since this is the periodic direction, 
whereas the dimensions along the x 2 and x 3 axes or the FG directions, h[ q \ h$\ and l[ r \ l 2 \ can 
vary from unit cell to unit cell. The dimensions of the subcells within a given cell along the FG 
directions are designated with running indices q and r which identify the cell number in the 
x 2 — *3 plane, where q and r remain constant along the x\ axis. For the remaining direction, x \ , 
the corresponding index p is introduced. Thus a given cell is designated by the triplet (p,q,r) for 
an infinite range of p due to periodicity in the x\ direction, and for q = 1,2, ..., N q and r = 1, 2, ..., 
N r , where N q and N r are the number of cells in the FG x 2 and x 3 directions. 

It is important to note that the unit cell (apy) in the present framework is not taken to be an 
RVE whose effective properties can be obtained through local homogenization, as is done in the 
standard uncoupled micromechanical approaches based on the concept of local action (Malvern, 
1969). In fact, for fully nonuniform distributions of fibers or inclusions in the x 2 -x 3 plane, no 
RVE can be identified. Thus the principle of local action is not applicable at the individual cell 
level, requiring the response of each cell to be explicitly coupled to the response of the entire 
array of cells in the FG directions. This is what is meant by the statement that the present 
approach explicitly couples the microstructural details with the global analysis, and thus sets 
HOTFGM-2D apart from the standard approaches found in the literature. The limitations of the 
standard uncoupled approach, and the error that results from decoupling of the local and global 
effects were recently discussed by Pindera et al. (1994, 1995). 
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2.1 Outline of the Solution Technique 

The solution of the thermo-mechanical boundary- value problem outlined in the foregoing is 
solved in two steps, following the general framework for the solution of the corresponding one- 
dimensional thermo-elastic problem discussed previously (Aboudi et al., 1993). In the first step, 
the temperature distribution in every cell is determined by solving the heat equation under 
steady-state conditions in each cell subject to the appropriate continuity and compatibility condi- 
tions. The solution to the heat equation is obtained by approximating the temperature field in 
each subcell of a unit cell using a quadratic expansion in the local coordinates x , x , x cen- 
tered at the subcell’s mid-point. A higher-order representation of the temperature field is neces- 
sary in order to capture the local effects created by the thermomechanical field gradients, the 
microstructure of the composite and the finite dimensions in the FG directions, in contrast with 
previous treatments involving fully periodic composite media which employed linear expansions 
(Aboudi, 1991). The unknown coefficients associated with each term in the expansion are then 
obtained by constructing a system of equations that satisfies the requirements of a standard 
boundary-value problem for the given temperature field approximation. That is, the heat equa- 
tion is satisfied in a volumetric sense, and the thermal and heat flux continuity conditions within 
a given cell, as well as between a given cell and adjacent cells, are imposed in an average sense. 

Given the temperature distribution in the functionally graded composite in the periodic and 
FG directions, internal displacements, strains and stresses are subsequently generated by solving 
the equilibrium equations in each cell subject to appropriate continuity and boundary conditions. 
The solution is obtained by approximating the displacement field in the FG directions in each 
subcell using a quadratic expansion in local coordinates within the subcell. The displacement 
field in the periodic x, direction, on the other hand, is approximated using linear expansion in 
local coordinates to reflect the periodic character of the composite’s microstructure along the x, 
axis. The unknown coefficients associated with each term in the expansion are obtained by satis- 
fying the appropriate field equations in a volumetric sense, together with the boundary condi- 
tions and continuity of displacements and tractions between individual subcells of a given cell, 
and between adjacent cells. The continuity conditions are imposed in an average sense. This 
results in a coupled system of equations involving the unknown coefficients in the displacement 
representation for each cell. 

An outline of the governing equations for the temperature and displacement fields in the 
individual subcells within the rows and column of cells considered in solving the outlined 
boundary-value problem is given in the following. A detailed derivation of these equations is 
presented in the Appendices A and B so as not to obscure the basic concepts by the involved 
algebraic manipulations. 
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2.2 Thermal Analysis: Problem Formulation 

Suppose that the composite material occupies the region lo^l < °°, 0 <x 2 ^H, 0 <x 3 <L. 
Let the composite be subjected to the temperature T t on the top surface (x 2 = 0), T B on the bot- 
tom surface (x 2 = H), T L on the left surface (x 3 = 0), and T R on the right surface (x 3 = L). Also, let 

N, 

N denote the number of cells in the interval 0 <x 2 <H, i.e., N q = H / £ (h\ q} + h%’). Likewise, 

q=\ 

N r 

let N r denote the number of cells in the interval 0 < x 3 < L, i.e., N r =L / £ (l\ r ^ + ). For q = 2, ... 

r = 1 

, N q _ i and r- 2, ... , N r -i the cells are internal, whereas for q = 1 ,N q and r = l, N r they are boun- 
dary cells. 

2.2.1 Heat Conduction Equation 

For a steady-state situation, the heat flux field in the material occupying the subcell (apy) of 
the (p,q,r)th cell, in the region defined by I -y^a> '*2 P> 1 '*3 1 - y4 r) ’ must 

satisfy: 

3 , q\ m + +hq$**=0 ( : 1 > 

where d l = d/dxf\ d 2 = d/dxf, d 3 = d/dxf. The components of the heat flux vector in this 
subcell are derived from the temperature field according to: 


= -k^diT^ , (i = 1, 2, 3; no sum) 


( 2 ) 


where J are the coefficients of heat conductivity of the material in the subcell (apy), and no 
summation is implied by repeated Greek letters in the above and henceforth. Given the relation 
between the heat flux and temperature, a temperature distribution that satisfies the heat conduc- 
tion equation is sought subject to the continuity and boundary conditions given below. 

2.2.2 Heat Flux Continuity Conditions 

The continuity of the heat flux vector q (aW at the interfaces separating adjacent subcells 
within the unit cell (p,q,r) is fulfilled by imposing the relations 


q\ l H> I 


ip.q.r) 

x?=d x /l 


-«f w 


, iP^r) 

1 


= -d 2 /l 


(3a) 
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(3b) 


i ^ r) i (p ' 9>r) 

«?'* I % .Ufa = <&* Ijf.-J*** 

(3c) 

In addition to the above continuity conditions within the (p,q,r)th cell, the heat flux continuity at 
the interfaces between neighboring cells must be ensured. The conditions that ensure this are 
given by 


. (p+Uq.r) . (p,q,r) 

q\ m 


(4a) 


(p,q+l,r) . (p.q.r ) 

<&*» k’ =*?>/* 


(4b) 


„ , (p.q.r+ 1 ) 


= <?^ 2) I *?=&/! 


(4c) 


2.2.3 Thermal Continuity Conditions 

The thermal continuity conditions at the interfaces separating adjacent subcells within the 
representative cell (p,q,r ) are given by relations similar to the corresponding heat flux continuity 
conditions. 


jKlpT) | 


(p.q.r) | 

jo , _ a t(2Pt) I -m . _ 

jc, =<i,/2 •* 1 JCi =-a 2 /2 


(P.q.r) 


(5a) 


l <P’9’ r > _ A | (P-9.'-) 

T (o1 ^ U'>=^/2 = ^ ls“=-*?>V2 


(5b) 


, (p.q.r) „ , (P.fl.r) 

T (aPI) I = I'pn = T (ap2) I — ^>/2 


(5C) 


while the thermal continuity at the interfaces between neighboring cells is ensured, as in the case 
of the heat flux field, by requiring that 
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(6a) 


. (p+\,q y r) 

T (m | M m _j i/2 = T {m 


. (P.9.'') 

I _ 0 > 

X] — ^2^ 


. (p,?+l,r) . (p,q,r) 

T {ali > I = r (a2lf) I g-wa ( 6b ) 

, (p y q,r+ 1) . (p,q,r ) 

r ( aPl) I ^ = _ /r , )/7 = t**® 1^=^ (6c) 


2.2.4 Boundary Conditions 

The final set of conditions that the solution for the temperature field must satisfy are the 
boundary conditions at the top and bottom, and left and right surfaces. The temperature in the 
cell (p,l,r) at the top surface must equal the applied temperature 7>, whereas in the cell (p,N q ,r ) 
at the bottom surface the temperature must be T B . 


I bri (1) 1 /iv 

rtoir>| =T t (x 3 ), x 2 =-±h[" (7a) 

. <P.N f .r) (2 ) 1 (N 

T {a2y) I =T b (jc 3 ), x 2 =jh^ (7b) 


where r = 1, • • ■ ,N r . 

Similarly, the temperature in the cell (p,q, 1) at the left surface must equal the applied tem- 
perature T l , whereas in the cell (p,q,N r ) at the right surface the temperature must be . 


y(api) 


(p.q. 1) 


= T l (x 2 ) , 




(8a) 


. (P.9.W 

j-(a$2) | 


Tr(x 2 ), 


-(2) 1 ,W, 


(8b) 


where q = 1, • • • ,N q . 

Alternatively, it is possible to impose mixed-boundary conditions involving temperature 
and heat flux at different portions of the boundary. 

2.3 Thermal Analysis: Solution 

The temperature distribution in the subcell (afty) of the (p,q,r)th cell, measured with respect 
to a reference temperature T re p is denoted by 7 ^°^. W e approximate this temperature field by a 
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second order expansion in the local coordinates xf*, xf\ and x^ as follows: 

T ( °M = + 3? TW + + y(3^i a)2 - ^-)TW> + {(3 xT - ^-)Tm 

1 /< r > 2 

+ |(33# )2 - J ^-)7$$ (9) 

where 7$8$, which is the temperature at the center of the subcell, and (l, m, n = 0, 1, or 2 
with / + m + n < 2) are unknown coefficients which are determined from conditions that will be 
outlined subsequently. It should be noted that eqn (9) does not contain a linear term in the local 
coordinates x^ . This follows directly from the assumed periodicity in the *i direction and sym- 
metry with respect to the lines xf* = 0 for a = 1 and 2. 

Given the six unknown quantities associated with each subcell (i.e., 7$$$, ... , 7$j§$) and 
eight subcells within each unit cell, 4SN q N r unknown quantities must be determined for a compo- 
site with N q rows and N r columns of different materials. These quantities are determined by first 
satisfying the heat conduction equation, as well as the first and second moment of this equation 
in each subcell in a volumetric sense in view of the temperature field approximation given by 
eqn (9). Subsequently, continuity of heat flux and temperature is imposed in an average sense at 
the interfaces separating adjacent subcells, as well as neighboring cells. Fulfillment of these field 
equations and continuity conditions, together with the imposed thermal boundary conditions at 
the top and bottom, and left and right surfaces of the composite, provides the necessary 4SN q N r 
equations for the 4SN q N r unknown coefficients in the temperature field expansion. We begin the 
outline of steps to generate the required 48 N q N r equations by first considering an arbitrary 
( p,q,r)th cell in the interior of the composite (i.e., q =2, ... , N q -l and r = 2, ... , A^-l). This pro- 
duces 4S(N q -2)(N r -2) equations. The additional equations are obtained by considering the boun- 
dary cells (i.e., q = 1, N q and r = 1, N r ). For these cells, most of the preceding relations also hold, 
with the exception of some of the interfacial continuity conditions between adjacent cells which 
are replaced by the specified boundary conditions. 

2.3.1 Heat Conduction Equations 

In the course of satisfying the steady-state heat equation in a volumetric sense, it is con- 
venient to define the following flux quantities: 
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( 10 ) 


QSfflL) = 


dji hfn. t?n 

J J f (xT )'(*? ) m 0f? Tq^d^cRfcEf 

-dji -hfn. -fin 


where l, m, n = 0 , 1 , or 2 with l +m+n< 2, and v[gf$ = d a h^ l\ r) is the volume of the subcell 
(apy) in the (p,q,r)th cell. For / = m =n =0, <2i($tf,o) is the average value of the heat flux com- 
ponent in the subcell, whereas for other values of ( l,m,n ) eqn (10) defines higher-order heat 
fluxes. These flux quantities can be evaluated explicitly in terms of the coefficients 7® by per- 
forming the required volume integration using eqns (2) and (9) in (10). This yields the following 
non-vanishing zeroth-order and first-order heat fluxes in terms of the unknown coefficients in 
the temperature field expansion: 


estito, 

(ID 


(12) 

eSSf.o, 

(13) 

eSffio) 

(14) 

(M2 

(15) 


Satisfaction of the zeroth, first and second moment of the steady-state heat equation results 
in the following eight relationships among the first-order heat fluxes i n the different sub- 

cells (aPy) of the (p,q,r)th cell, after some involved algebraic manipulations (see Appendix A): 

t eMfojAfc + oM + xA =o (i6) 


where the triplet (apy) assumes all permutations of the integers 1 and 2. 

2.3.2 Heat Flux Continuity Equations 

The continuity of heat fluxes at the subcell interfaces associated with the periodic x x direc- 
tion, eqn (3a) imposed in an average sense, is ensured by: 
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[ ] <p,q ' r) = 0 


(17) 


We note that eqn (4a) that ensures continuity of heat flux in the x\ direction between neighbor- 
ing cells is identically satisfied by the chosen temperature field representation due to the periodi- 
city of the composite in this direction. 

The equations that ensure heat flux continuity at the subcell interfaces, as well as between 
individual cells, associated with the x 2 and x 3 directions, eqns (3b), (3c) and eqns (4b), (4c) are 
given by 


[ -12 Q%$, 0) /h , + <$X) - <>Q$b)/h2 ]** r) - [ figft 0) + = 0 (18) 

[ ~Q$l to) + \( Q%to) - 6Q^lo)/h 2 )] ( ™' ) + yl Q&o) + 6<2^,0)/fc 2 ] (p ’*' , ’ r) =° (19) 

[ -120^,,)/Zi + Qm,0) ~ 60$ 2 o!i)/?2 ] (p ’ 9 ’ r) - 1 Q$$o) + 6Q$ 2 J 1) /l 2 l <p ’* r - U = 0 (20) 

[ -fig® 0) + {( 69.0) - 60 ^i>// 2 )] (p ’ 9 ’ r) + yt 0§&) + 6fig® i^ 2 = 0 ( 21 ) 


Equations (17) - (21) provide us with twenty additional relations among the zeroth-order 
and first-order heat fluxes. These relations, together with eqn (16), can be expressed in terms of 
the unknown coefficients by making use of eqns (11) - (15), providing a total of twenty- 
eight of the required forty-eight equations necessary for the determination of these coefficients 
in the (p,q,r)th cell. 

2.3.3 Thermal Continuity Equations 

An additional set of twenty equations necessary to determine the unknown coefficients in 
the temperature field expansion is subsequently generated by the thermal continuity conditions 
imposed on an average basis at each subcell and cell interface. Imposing the thermal continuity 
at each subcell interface in the periodic x t direction, eqn (5a), we obtain the following conditions 
for the (p,q,r)th cell: 


[ T\U + \d\tM t' q ' r) = [ T?M + \d\T%U ^' q ' T) 


( 22 ) 


We note that the continuity of temperature between neighboring cells in the x x direction, eqn 
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(6a), is automatically satisfied by the chosen temperature field representation which reflects the 
periodic character of the composite in this direction. 

The continuity of temperature at the interfaces between the subcells, as well as between the 
neighboring cells, in the FG directions, eqns (5b) - (5c) and eqns (6b) - (6c), on the other hand, 
yield 



Equations (22) - (26) comprise the required additional twenty relations. 

2.3.4 Governing Equations for the Unknown Coefficients in the Temperature Expansion 

The steady-state state heat equations, eqn (16), together with the heat flux and thermal con- 
tinuity equations, eqns (17) - (21) and eqns (22) - (26), respectively, form altogether 48 linear 
algebraic equations which govern the 48 field variables in the eight subcells (ccpy) an 
interior cell (p,q,r ); q = 2, ... , N q - 1, r = 2, ... , N r — 1. For the boundary cells q = lj'fq, and r = 1 7V r , 
a different treatment must be applied. For q = 1, the governing equations, eqns (16) - (17) and 
(20) - (26) are operative. Relations (18) - (19), on the other hand, which follow from the con- 
tinuity of heat flux between a given cell and the preceding one are not applicable. They are 
replaced by the condition that the heat flux at the interface between subcell (aly) and (oc2y) of the 
cell (p, l,r) is continuous, as well as the applied temperature relation at the surface x 2 = 0, eqn 
(7a). For the cell q = N q , the previous equations are applicable except eqns (24) which are obvi- 
ously not operative. These equations are replaced by the specific temperature applied at the sur- 
face x 2 = H, eqn (7b). Similar reasoning holds for the subcells r = 1 and r = N r . 

The governing equations at the interior and boundary cells form a system of 48 N q N r linear 
algebraic equations in the unknown coefficients . Their solution determines the temperature 
distribution within the FG composite subjected to the boundary conditions (7) and (8). The final 
form of this system of equations is symbolically represented below 
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kT = t 


( 27 ) 


where the structural thermal conductivity matrix k contains information on the geometry and 
thermal conductivities of the individual subcells (apy) in the N q N r cells spanning the x 2 and x 3 
FG directions, the thermal coefficient vector 7 contains the unknown coefficients that describe 
the thermal field in each subcell, i.e., 7 = [ 7$$, .... , ] where 7$™$ = [ 7(000), f( 010 ). 7(001), 

T( 200 ), 7 ( 020 ) 5 7 ( 002 ) ] (otP7) , and the thermal force vector t contains information on the boundary 
conditions. 

2.4 Mechanical Analysis: Problem Formulation 

Given the temperature field generated by the applied temperatures T t , Tg, and T L , T R 
obtained in the preceding section, we proceed to determine the resulting displacement and stress 
fields. This is carried out for arbitrary mechanical loading applied to the surfaces of the compo- 
site in the ;c 2 -*3 plane, excluding shearing in the x, direction. 

2.4.1 Equations of Equilibrium 

The stress field in the subcell (aPy) of the (p,q,r)th cell generated by the given temperature 
field must satisfy the equilibrium equations 

a, off* + a 2 o^> + a 3 a£ w = o , j = 1, 2, 3 (28) 


where the operator 9, has been defined previously. The components of the stress tensor, assum- 
ing that the material occupying the subcell (aPy) of the (p,q,r)th cell is orthotropic, are related to 
the strain components through the familiar generalized Hooke’s law: 

- r^7 (aW (29) 

where are the elements of the stiffness tensor and the elements of the so-called ther- 
mal tensor are the products of the stiffness tensor and the thermal expansion coefficients. The 
components of the strain tensor in the individual subcells are, in turn, obtained from the strain- 
displacement relations 

ej = y (9 + djtfi**) , i,j = 1 , 2, 3 (30) 
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Given the relation between the stresses and displacement gradients obtained from eqns (29) and 
(30), a displacement field is sought that satisfies the three equilibrium equations together with 
the continuity and boundary conditions that follow. 

2.4.2 Traction Continuity Conditions 

The continuity of tractions at the interfaces separating adjacent subcells within the repeat- 
ing unit cell (p,q,r) is fulfilled by requiring that 


■ (p,q,r) . {p>q,r ) 

(31a) 

( p,q,r ) . (p,q,r) 

(31b) 

, (p&r) , (p.q,r) 

off" \^= ir/2 =off 2) Ur = -/-/2 

(31c) 


In addition to the above continuity conditions within the (p,q,r)th cell, the traction continuity at 
the interfaces between neighboring cells must be ensured. These conditions are fulfilled by 


requiring that 


. (p+hq,r) n i (p,qs) 

1 tf=- dl/ 2=o[?» 1 

(32a) 

(p,?+l,r) , (P.q.r) 

off* \^^n=off* 1 

(32b) 

, {p,q,r+\) , (p,q>r) 

off" 1 $=-w*=off 2) \g.tpn 

(32c) 


2.4.3 Displacement Continuity Conditions 

At the interfaces of the subcells within the repeating unit cell {p,q,r ) the displacements 
u = (u], h 2 , w 3 ) must be continuous, 


u m | 

u wy) | 

„(opl) I 


(p,q.r) 

^ 1> =d 1 /2=« 


(p,q.r) 


(p,q,r) 

x" = /y >/2 = “ 


(2pTf) | 
(o2-y) | 
(«P2) | 


(p,q,r) 

_» . . 

•*1 = -difl 


(P,q,r) 

sf =-/>?>/ 2 


(p.q.r) 

X? =-i¥n 


(33a) 

(33b) 

(33c) 
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while the continuity of displacements between neighboring cells is ensured by requiring that 


u 

H (allO 

M (cepi) 


(p+l,q,r) , (p,q,r ) 

_a> . „ = I -jS> 

a:, = — u ' a:, 


<f,/2 


(P.9+1. ') . _ . i (P.?.' - ) 

- = u (a2i) I 






(p.?,r+l) 

x" =-V /l 


„ , (P.f.r) 

— U («P» I _W 

“ * *3 = 




(34a) 

(34b) 

(34c) 


2.4.4 Boundary Conditions 

The final set of conditions that the solution for the displacement field must satisfy are the 
boundary conditions at the top and bottom, and left and right surfaces. The traction vector in the 
cells (p, l,r) and (p,N q ,r ) at the top and bottom surfaces, respectively, must equal the applied sur- 
face loads, 


(p. l,r) 

I =/ri(*3)> 


( P.N,,r ) 

off* I = /*(* 3) , 


4” 


(35a) 

(35b) 


where r = f Ti and f Bi describe the spatial variation of these loads at the top and bottom 

surfaces. Similarly, if the right or left surfaces are rigidly clamped (say), then 


„(aPD | 
M (“P 2 ) | 


( P.9 . 1) 

= 0, 


(P.9.W r) 

= 0, 


4 2) = }4 W 


(36a) 

(36b) 


where q = 

For other types of boundary conditions, eqns (35) - (36) should be modified accordingly. 

2.5 Mechanical Analysis: Solution 

Due to symmetry considerations, the displacement field in the subcell (afty) of the (p,q,r)th 
cell is approximated by a second-order expansion in the local coordinates * 1 °°, xf \ and as 
follows: 
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uf *» -jfT’wlsBfe 


= wga> +^ p> ^«) +*?w«> + {<3*T 2 - + 

i(3*?® - |*jf l> )WSSj) + }(3I?’ 2 - ■j-ir)"*) (37) 

«?W = Wgffi) +*? wSffi, +I?W®1, + |(3lf“ - j4)W$ffi> + 

|(3lf 2 - )»$$&> + y(3J? 2 - 


where Wj“<§$), which are the displacements at the center of the subcell, and W$J&) (i = 1,2,3) 
must be determined from conditions similar to those employed in the thermal problem. In this 
case, there are 104 unknown quantities in a subcell (apy). The determination of these quantities 
parallels that of the thermal problem. Here, the heat conduction equation is replaced by the three 
equilibrium equations, and the continuity of tractions and displacements at the various interfaces 
replaces the continuity of heat fluxes and temperature. Finally, the boundary conditions involve 
the appropriate mechanical quantities. As in the thermal problem, we start with the internal cells 
and subsequently modify the governing equations to accommodate the boundary cells q — \, N q , 
and r = 1, N r . 

It should be noted that the first equation in (37) does not contain linear terms in the local 
coordinates x? ] and xf . This follows from the assumed periodicity in the xj direction and sym- 
metry with respect to x\ — 0 (cc = 1 , 2 ). Further, the absence of a constant term in the first equa- 
tion that represents subcell center x \ - displacements, say leads to the result that the aver- 

age normal strain of the composite associated with the x i direction is zero. It is possible to gen- 
eralize the present theory by including subcell center displacements associated with the * i direc- 
tion that produce uniform composite strain i . This generalization leads to an overall behavior 
of a composite, functionally graded in the X 2 and *3 directions, that can be described as a gen- 
eralized plane strain in the x\ direction. This generalization is not trivial as it requires coupling 
between the present higher-order theory and an RVE-based theory which employs a homogeni- 
zation scheme (Aboudi, et al„ 1995). In Section 2.5.5 we briefly outline how the present formu- 
lation can be modified to admit generalized plane strain in the x x direction in order to be able to 
carry out comparison between HOTFGM-2D and finite-element analysis. 
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2.5.1 Equations of Equilibrium 

In the course of satisfying the equilibrium equations in a volumetric sense, it is convenient 
to define the following stress quantities: 



For / = m = n = 0, eqn (38) provides average stresses in the subcell, whereas for other values of 
(l,m,n) higher-order stresses are obtained that are needed to describe the governing field equa- 
tions of the higher-order continuum. These stress quantities can be evaluated explicitly in terms 
of the unknown coefficients Wjflgg) by performing the required volume integration using eqns 
(29), (30) and (37) in eqn (38). This yields the following non-vanishing zeroth-order and first- 
order stress components in terms of the unknown coefficients in the displacement field expan- 
sion: 

sl*„) (39) 

sMi,o> = 8 , - (40) 

sf#o.i> = (41) 


with similar expressions for S^^o,o>* S^ftlYo)- SSfSVi). and ^n( 8 ?o, 0 ). ^CrfcCi.Oji ^ 3 ^?o.d> and 



(42) 


(43) 

5^0,0) =cSW(WS®l) + 

(44) 


(45) 

s&fto. D = |4 r)2 c^w^) 

(46) 
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Satisfaction of the equilibrium equations results in the following sixteen relations among 
the volume-averaged first-order stresses in the different subcells (<x(3y) of the 

cell, after lengthy algebraic manipulations (see Appendix B): 

[ »/<£ + sgw mH + = 0 < 47 > 

where j = 2, 3 and, as in the case of eqn (16), the triplet (afty) assumes all permutations of the 
integers 1 and 2. 

2.5.2 Traction Continuity Equations 

The continuity of tractions at the subcell interfaces associated with the periodic *1 direc- 
tion, eqn (31a) imposed in an average sense, is ensured by the following relations: 

[ 5*0.0) -S®o.o>] ( ' w >=0 (48) 

[ i + S®W«2 = 0 (49) 

[ S®lW«f , + ] <w> =° < 50 > 

We note that eqn (32a), which ensures continuity of tractions between adjacent cells in the 
periodic x x direction, is identically satisfied by the chosen displacement field representation due 
to the periodic character of the composite material in this direction. 

The equations that ensure traction continuity between individual subcells, as well as 
between individual cells, associated with the *2 3nd x$ directions, eqns (31b), (31c), and eqns 
((32b), (32c), are given by 


[-IMg'Jtl.oi'fci +SgSfo.o) 

[ -Sffim + 0.0) - 3S$*,.o,/». 2 l g> ' , ' ,) + {l 5».0) 

[-125gtVi)f>i +SSS®U 
[ -Sgffilxo, + |s$S).0> - 3Sgg'o.l)/>2 4 {[ 


+6igafi.o)/h2i ip - , -'' ,, =o 
+ 6Sgo!fi,0)/J'2l < '' , -‘-'>=0 

+ 6SgS5J,,„/l2 ]<w-»=0 


(51) 

(52) 

(53) 

(54) 


where j = 2 and 3. 
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Equations (48) - (54) provide us with forty-four (44) additional relations among the zeroth- 
order and first-order stresses. These relations, together with eqn (47), can be expressed in terms 
of the unknown coefficients by making use of eqns (39) - (46), providing a total of sixty 

(60) of the required one hundred and four (104) equations necessary for the determination of 
these coefficients in the (p,q,r)th cell. 

2.5.3 Displacement Continuity Equations 

The additional forty-four (44) relations necessary to determine the unknown coefficients in 
the displacement field expansion are subsequently obtained by imposing displacement continuity 
conditions on an average basis at each subcell and cell interface. The continuity of displacements 
in the periodic x x direction at each subcell interface of the (p,q,r)th cell, eqn (33a), is satisfied by 
the following conditions: 


+ d 2 W[ 2 <$ ) ]<*«•'> =0 (55) 

[ ~ ]<**'> = 0 (56) 

[ - Wf(68» - jdM 2 ^ ]<*«•'> = 0 (57) 


The displacement continuity between neighboring cells in the periodic x x direction, eqn (34a), is 
automatically satisfied by the chosen displacement field representation which reflects the 
periodic character of the composite in this direction. 

The displacement continuity conditions at the inner interfaces, eqns (33b) - (33c), as well as 
between neighboring cells, eqns (34b) - (34c), in the FG x 2 and x 3 directions, yield 
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Equations (55) - (65) comprise the required additional forty-four relations. 

2.5.4 Governing Equations for the Unknown Coefficients in the Displacement Expansion 

The equilibrium equations, eqn (47), together with the traction and displacement continuity 
equations, eqns (48) - (54) and eqns (55) - (65), respectively, form altogether 104 equations in 
the 104 unknowns which govern the equilibrium of a subcell (a(3y) within the (p,q,r)th 

cell in the interior. As in the thermal problem, a different treatment must be adopted for the 
boundary cells ip , l,r), (p,N q ,r), and (p,q, 1) and (p,q,N r ). For (p, l,r), the above relations are 
operative, except eqns (51) and (52), which follow from the continuity of tractions between a 
given cell and the preceeding one. These eight equations must be replaced by the conditions of 
continuity of tractions at the interior interfaces of the cell (p, l,r) and by the applied tractions at 
*i =0, eqn (35a). For the cell (p,N q ,r ), the previously derived governing equations are operative 
except for the four relations given by eqns (59) and (61) which are obviously not applicable. 
These are replaced by the imposed traction conditions at the surface x 2 = H, eqn (35b). Similar 
arguments hold for boundary cells (p,q, 1) and (p,q,N r) . Consequently, the governing equations at 
both interior and boundary cells form a system of \04N q N r linear algebraic equations in the field 
variables within of the cells of the functionally graded composite. The final form of this system 
of equations is symbolically represented below 

KU=f (66) 

where the structural stiffness matrix K contains information on the geometry and thermo- 
mechanical properties of the individual subcells (a|3y) within the cells comprising the function- 
ally graded composite, the displacement coefficient vector U contains the unknown coefficients 
that describe the displacement field in each subcell, i.e., U = [ .... , U { ] where = 
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[ Wi(ioo), • • • > ^ 3 ( 002 ) ] (0 ^ Y) . die mechanical force vector /contains information on the boun- 
dary conditions and the thermal loading effects generated by the applied temperature. 

2.5.5 Extension to Generalized Plane Strain in the xj Direction 

As indicated previously, the present formulation in which the displacement components 
within the subcell (a(5y) of the (p,q,r)th cell are expanded in accordance with eqn (37) leads to 
the result that the strain averaged over the entire volume V occupied by the functionally 
graded composite is zero. 


X i «^ = 0 (67) 

V q= 1 r = 1 a,p, 1 *=l 

It is frequently desirable to obtain a generalized plane strain situation where e u = non-zero con- 
stant. This can be achieved by adding a constant term to the first equation in (37) and applying a 
homogenization procedure in the x x periodic direction. It can be shown that the only equation 
affected by this process is eqn (55) which must be replaced by (Aboudi, et al., 1995): 


[dM l M» +d 2 W$b ]<™ r) = (<*i +d 2 )e» 


( 68 ) 


where Eu is an unknown value that must be determined from the condition that On - 0, where 
On is the average of over the entire volume of the composite: 

5i,=vX X £ vggoftW-O (69) 

y q = l r= 1 


3.0 APPLICATIONS 

The approach outlined in the foregoing is employed to investigate the response of a sym- 
metrically laminated B/Ep-Ti plate subjected to a uniform temperature change of -154.45°C. 
This temperature change simulates cool down from the fabrication temperature which induces 
residual stresses into the individual plies due to a thermal expansion mismatch between the 
boron/epoxy and the titanium plies. The residual stress fields exhibit large gradients near the free 
edges of the laminate which were investigated by Herakovich (1976) using the finite-element 
approach. Herein, we first compare the finite-element results for the stress fields near the free 
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edge with the predictions of HOTFGM-2D, treating the boron/epoxy (B/Ep) plies as homogene- 
ous with equivalent effective (or homogenized) properties. This comparison demonstrates the 
capability of the new coupled theory to capture large gradients in the stress fields at geometric 
and material discontinuities (i.e., along interfaces at the free edge of a laminated plate) in the 
presence of a spatially uniform temperature field. Subsequently, the effect of microstructure of 
the B/Ep plies on the free-edge stress fields is investigated in view of the large fiber diameter of 
the boron fibers and relatively small thickness of the B/Ep plies. Finally, the utility of nonuni- 
formly distributing (i.e., functionally grading) boron fibers in the B/Ep plies near the free edge to 
reduce the large stress gradients in this region is demonstrated. 

The cross-section geometry of the investigated laminate is given in Figure 3. The thickness, 
designated by H in the figure, and the width L produce a laminate with an aspect ratio of L/H = 
12.5. The direction of the boron fibers in the external B/Ep plies is parallel to the (out-of- 
plane) axis, along which the laminate is considered to be infinitely long. The volume fraction of 
the fibers is 0.50. The resulting macroscopic thermo-mechanical properties of the B/Ep plies and 
the titanium inner sheets are given in Table 1. As in Herakovich (1976), these properties are con- 
sidered to be temperature-independent. Comparison of the thermal expansion coefficients of the 
B/Ep plies and the titanium sheets reveals a significant mismatch in the transverse direction 
(along the x 3 axis), with a smaller mismatch in the out-of-plane direction. Figure 4 illustrates the 
axial and transverse thermal expansion coefficients, and (Xf, respectively, of the B/Ep ply as a 
function of the fiber content Vf generated using the method of cells micromechanics model 
(Aboudi, 1991). Included in the figure is the thermal expansion coefficient of the isotropic 
titanium sheet. The thermoelastic properties of the boron fibers and the epoxy matrix used to 
generate this figure are given in Table 2. The graphical results shown in Figure 4 indicate that 
the large thermal expansion mismatch between the B/Ep and the Ti plies in the transverse direc- 
tion (i.e., jc 3 ) can be reduced by increasing the fiber content of the B/Ep ply above the current 
value of 0.50. It is expected that the reduction of the transverse CTE mismatch will lead to a 
decrease in the transverse residual stresses and thus the severity of the free-edge stress gradients. 
Figure 4 thus provides a motivation for using functionally graded fiber architectures in the vicin- 
ity of the free edge to decrease the transverse thermal expansion mismatch in this region, and 
thus reduce the interlaminar stresses. Since the interlaminar stress field is a localized effect, it is 
reasonable to expect that it is sufficient to limit grading of the fiber content in the B/Ep plies to 
the immediate vicinity of the free edge. 
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3.1 Comparison of HOTFGM-2D and FE Predictions Based on Homogenized B/Ep Pro- 
perties 

Figure 5 presents comparison of the normal (on » ri 22 > O 33 ) *®d shear (o 23 ) stress distribu- 
tions in the titanium sheet along the interface separating the B/Ep and Ti plies (see Figure 3) 
generated using HOTFGM-2D and the finite-element analysis of Herakovich (1976). The 
fini te-element results were generated using three different meshes, indicated by B-l, B-2 and B- 
3 in the figure, with each successive mesh undergoing increasingly greater refinement. The 
HOTFGM-2D and finite-element results were obtained using the effective or homogenized 
thermo-mechanical properties of the B/Ep and Ti plies given in Table 1. As stated previously, 
the illustrated stress distributions were induced by subjecting the [B/Ep-Ti]* laminate to a uni- 
form temperature change of -154.45°C. Due to the symmetry of the considered plate with respect 
to the x 2 and x 3 axes, only one quarter of the plate needs to be analyzed (see Figure 3) under 
appropriate boundary conditions which reflect these symmetries. 

As observed in the figure, the stress profiles exhibit large stress gradients in the immediate 
vicinity of the free edge. Away from the free edge, these stress distributions attain uniform 
values that can be predicted using the classical lamination theory for sufficiently large L/H 
aspect ratios. The rapid decay of the interlaminar stresses to their lamination theory or far-field 
values occurs over a distance that is approximately one laminate thickness H from the free edge. 
The far-field values of the normal stresses o n , o 22 and o 33 predicted by the coupled higher-order 
theory are 29.9, 0.0, and -54.95 MPa, respectively. The far-field value of the shear stress o 23 is 
0.0. These results coincide with the lamination theory predictions. 

In the vicinity of the free-edge, the large gradient and magnitude of the o n stress com- 
ponent obtained from the coupled higher-order analysis compares very favorably with the 
finite-element results generated with the most refined mesh. The behavior of the <r 22 stress com- 
ponent near the free edge predicted by the higher-order theory also compares favorably with the 
finite-element results, lying between the B-l and B-3 mesh predictions. This component is often 
responsible for delamination initiation at the free edge when it is tensile, as is the present situa- 
tion. The normal stress component a 33 is in the direction of the free edge and thus has to vanish 
on the lateral surface x 3 /L = 0.5. Both the higher-order theory and finite-element predictions 
indicate that this stress component does tend to zero with decreasing distance from the from 
edge. The finite-element results indicate an initial decrease in this stress component relative to 
the far-field value (i.e, an increase in the magnitude of the compressive stress) followed by rapid 
reversal and decay to zero. The magnitude of this initial decrease predicted by the HOTFGM- 
2D analysis however, is substantially smaller relative to the finite-element predictions. On the 
other hand, the magnitude of the normal stress o 33 in the immediate vicinity of the free edge is 
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closer to zero than that predicted by the finite-element analysis using the most refined mesh (B- 
3). Finally, the comparison of the shear stress components 023 predicted by the two approaches is 
also favorable. This stress component also undergoes a rapid reversal at the free edge, initially 
decreasing (i.e., increasing in magnitude), then reversing direction and rapidly decaying to zero. 
The magnitude of the maximum shear stress at the reversal point predicted by the coupled 
higer-order theory is somewhat smaller than that predicted by the finite-element analysis. On the 
other hand, the shear stress at the free edge predicted by HOTFGM-2D is much smaller (in fact 
almost zero) than that predicted by the finite-element approach. 

The comparison of the stresses in the titanium sheet at the B/Ep-Ti interface generated 
using the HOTFGM-2D and finite-element approaches indicates that the major features of the 
near free-edge stress fields are correctly captured by the coupled higher-order theory. While in 
some instances the actual magnitudes are not in perfect agreement, it is not clear at this point 
whether the problem lies with the HOTFGM-2D or finite-element results since the discrepan- 
cies, when they occur, are not sufficiently consistent to point to either of the two approaches as a 
potential culprit. It is clear however, based on the presented comparison, that the coupled 
higher-order theory is sufficiently sensitive to capture the large stress gradients, and even rapid 
stress reversals (see <j 23 distribution in Figure 5) that occur in regions of geometric discontinui- 
ties such as the free edge. This sets the stage for investigating the effect of microstructure on the 
free-edge stress fields. 

3.2 Effect of Microstructure of the B/Ep Plies on the Free-Edge Stress Fields 

The results presented heretofore have been generated by treating the B/Ep plies as homo- 
geneous with certain effective or homogenized thermo-mechanical properties. These properties 
can be either measured directly in the laboratory or determined using a micromechanics scheme 
when macroscopically homogeneous deformation fields are imposed. However, when the rein- 
forcement size (i.e., fiber diameter) is large relative to the thickness of a composite ply, which is 
the case with boron-reinforced and silicon carbide-reinforced composites, the meaning of a 
material property becomes fuzzy in the presence of large stress gradients as previously dis- 
cussed. This is the case in the present situation at the free edge given the large boron fiber diam- 
eter relative to the thickness of the B/Ep ply. It is thus important to characterize the error intro- 
duced in the analysis of free-edge stress fields based on the homogenized properties of the B/Ep 
ply. 

Figure 6 illustrates eight configurations of the [B/Ep-Ti], laminate with increasingly refined 
microstructure of the B/Ep plies at the free edge that were investigated in order to determine the 
effect of the microstructure on the stress fields. The overall dimensions of the laminate and the 
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individual plies are the same as in the previous problem with homogenized B/Ep properties. The 
number of boron fibers in the through-thickness direction of the B/Ep plies was taken as two, 
and the fiber dimensions were chosen to yield the same fiber volume fraction of 0.50 as in the 
previous problem. The properties of the boron fiber and epoxy matrix reported in Table 2 pro- 
duce the same homogenized properties for the B/Ep plies as those employed previously. These 
properties were assigned to the fiber and matrix phases in the regions of the B/Ep plies shown in 
Figure 6 with increasingly refined microstructures. Outside of those regions, homogenized 
thermo-mechanical properties were employed for the B/Ep plies. 

Figure 7 presents the normal stress a 2 2 in the titanium sheet at the B/Ep-Ti interface for the 
eight configurations shown in Figure 6. This stress component, due to its tensile character, is 
responsible for the initiation of delamination, and thus is of particular importance for the con- 
sidered configuration subjected to the given thermal load. Examination of the individual figures 
for each configuration reveals an increasingly complex character of the stress distribution near 
the free edge with increasing refinement of the microstructure that is due to the interaction with 
the individual boron fibers in the adjacent B/Ep ply. The stress distributions are characterized by 
rapid oscillations that coincide with the locations of the boron fibers in the B/Ep plies directly 
above the titanium sheet. The reversals in the sign of the stress oscillations could potentially 
have an effect on the delamination process. More importantly, however, the maximum normal 
stress at the free edge increases with increasing refinement of the microstructure. This is clearly 
illustrated in Figure 8 where the maximum stress a 22 at the free edge for the eight different con- 
figurations has been normalized by the corresponding stress obtained from the HOTFGM-2D 
analysis using homogenized B/Ep properties. As is observed, the actual maximum stress at the 
free edge in the presence of microstructure asymptotically approaches a uniform value with 
increasing number of cells in the B/Ep ply at the free edge. At least seven cells are required in 
the horizontal direction at the free edge to capture the actual maximum value of in the 
titanium ply. More importantly, the actual maximum stress is 35% greater than the correspond- 
ing stress obtained using homogenized properties of the B/Ep plies. This is a significant differ- 
ence that cannot be overlooked in predicting the onset of delamination, revealing the shortcom- 
ing of the homogenized continuum approach in the presence of course microstructure and large 
stress gradients. 

3.3 Delamination Control Through Microstructural Tailoring of the B/Ep Plies 

As illustrated in the preceding sections, the large thermally-induced interlaminar stresses in 
the vicinity of the free edge, and in particular the large tensile normal g 22 stress commonly 
called peel stress , may be sufficiently large to initiate delamination during fabrication cool down 
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or subsequent mechanical loading. It is of technological interest to reduce this stress component 
in order to increase the load-bearing capability of such laminates. In this section we investigate 
the possibility of accomplishing this through the use of functionally graded architectures in the 
B/Ep plies. 

We choose two approaches to reduce the peel stress at the free edge. In the first approach, 
we selectively remove fibers from the B/Ep plies in the vicinity of the free edge in order to 
create a local clamping in the vertical direction at the lateral surface of the laminate. It is 
presumed that this localized clamping arises from the tendency of the matrix phase to contract 
more than the surrounding B/Ep material under the given temperature change, and that the con- 
traction in the vertical direction has a greater effect on the stress field than the contraction in the 
horizontal direction. Since eight fiber/matrix cells in the B/Ep plies at the free edge are sufficient 
to capture the actual stress field in the presence of microstructural details (see Figure 8), we use 
this configuration to test the hypothesis described above. Figure 9 presents the interlaminar peel 
stress profiles generated by removing one and two columns of fibers in the vicinity of the free 
edge. Figure 9a illustrates the interlaminar peel stress distribution for the baseline eight-cell con- 
figuration with no columns of fibers removed, while Figures 9b, 9c and 9d present the 
corresponding distributions for the configurations with column 2, columns 3 and 4, and columns 
2 and 3 removed. The results indicate that the removal of the columns of fibers generally tends 
to lower the normal peel stress in the regions directly below the missing fibers. Further, the mag- 
nitude of the compressive stress in the cells adjacent to the free edge is also increased. However, 
the maximum tensile stress at the free edge itself is actually increased when the fibers are 
removed. This is illustrated in Figure 10 which presents the maximum free-edge peel stress in 
the configurations with the missing fibers normalized with respect to the corresponding value of 
the baseline configuration given in Figure 9a. The greatest increase in the maximum value of the 
free-edge peel stress occurs in the configuration with the columns 2 and 3 removed, while the 
smallest increase occurs in the configuration with the columns 3 and 4 removed. Thus it appears 
that the increased tendency of the matrix with the missing fibers to contract in the vertical direc- 
tion is more than offset by the horizontal contraction tendency. The overall effect is to increase 
the thermal expansion mismatch between the B/Ep and Ti plies in the horizontal direction, 
resulting in the increase of the maximum peel stress at the free edge. 

The above exercise also sheds light on the effect of imperfections in fiber spacing near the 
free edge that may arise due to poor fiber placement control during fabrication. The results sug- 
gest that a localized increase in the thermal expansion mismatch between adjacent plies caused 
by missing fibers near the free edge will have a detrimental effect on the laminate’s delamina- 
tion resistance. 
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In the second approach used to reduce the interlaminar peel stress at the free edge, the fiber 
spacing in the horizontal direction was decreased in a linear manner with decreasing distance 
from the free edge. This effectively increases the local fiber volume fraction in the B/Ep plies 
which, in turn, decreases the transverse thermal expansion mismatch between the adjacent plies 
at the free edge as is observed in Figure 4. The results presented in Figures 9 and 10 indeed indi- 
cate that the local fiber volume fraction has a significant effect on the interlaminar stresses at the 
free edge, thereby supporting this second approach in reducing the peel stress. In addition to 
decreasing the horizontal fiber spacing in the vicinity of the free edge, the two rows of fibers 
were also shifted vertically in a uniform manner in order to bring them closer to the B/Ep-Ti 
interface and thus decrease the local thermal expansion mismatch in the thickness direction. 

The functional grading of the fiber distribution in the horizontal direction was accom- 
plished as follows. The total horizontal distance occupied by the eight fiber cells was 2032 pm. 
The distance of the center of the first fiber from free edge was taken to be 97.36 pm. This was 
chosen such that the distance between the eighth fiber and the homogenized B/Ep material was 
215.9 pm. This is exactly one half of the width of the unit cell with a fiber volume fraction of 
0.50 in the uniformly spaced configuration. The distances between the centers of the eight fibers 
were linearly increased with increasing distance from the free edge in the manner illustrated in 
Table 3. 

Figure 1 1 presents the four different fiber architectures near the free edge in the B/Ep plies 
that were generated using the combination of horizontal and vertical shifting of the boron fibers. 
The relative locations of the fibers with respect to the dimensions of the B/Ep plies are the same 
as in the actual configurations. The baseline configuration is the uniformly spaced configuration 
considered previously with center-to-center fiber spacing of 254.0 pm, shifted horizontally 
towards the free edge such that the location of the center of the first fiber coincides with the 
center of the first fiber in the linearly spaced configuration (i.e., 97.36 pm from the free edge). 
The second configuration was generated from the first by uniformly shifting the two rows of 
fibers in the vertical direction while preserving the uniform horizontal fiber spacing. The third 
configuration was generated by linearly increasing the fiber spacing in the horizontal direction 
with increasing distance from the free edge, in the manner described previously, while preserv- 
ing the vertical spacing of the baseline configuration. Finally, the fourth configuration was 
obtained from the third by vertically shifting the two rows of fibers in the same manner as was 
used to generate the second configuration. The total fiber volume fraction of the B/Ep region 
with the four fiber architectures considered was 0.50. The local fiber volume fractions of the 
vertical columns of cells in the uniformly spaced and linearly spaced fiber configurations are 
given in Table 4. 
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The resulting peel stress distributions along the B/Ep-Ti interface in the four configurations 
caused by the uniformly applied temperature change of -154.45°C are illustrated in Figure 12. 
Comparing the peel stress distribution in the uniformly-spaced baseline configuration, Figure 
12a, with the stress distribution in the second configuration. Figure 12b, we observe increased 
peel stress oscillations in the regions directly below the fiber locations. The relative increase in 
the oscillations is caused by the reduced distance between the interface and the bottom row of 
fibers produced by the vertical shift. The vertical shift also produces a reduction in the maximum 
peel stress at the free edge relative to the baseline configuration. The peel stress distribution in 
the configuration with the linearly spaced fibers, Figure 12c, exhibits decreasing oscillations 
with decreasing distance from the free edge relative to the baseline configuration. A substan- 
tially greater decrease in the maximum peel stress at the free edge relative to the baseline confi- 
guration is also observed, as compared to the free-edge peel stress reduction observed in the 
second configuration. Finally, the peel stress distribution in the fourth configurations presented 
in Figure 12d exhibits characteristics that are peculiar to the second and third configurations. 
Here, we observe an increase in the peel stress oscillations relative to the baseline configuration 
that decrease with decreasing distance from the free edge. The reduction in the maximum peel 
stress at the free edge also appears to be substantial, and slightly greater than in the third confi- 
guration. A more precise comparison of the maximum peel stress values at the free edge in the 
functionally graded configurations is presented in bar chart format in Figure 13, normalized by 
the corresponding value of the uniformly-spaced baseline configuration. As is observed, the 
greatest reduction in the peel stress (i.e., approximately 25%) occurs when the fibers are both 
linearly spaced and vertically shifted. When the fibers are linearly spaced without the vertical 
shift, the reduction is approximately 24%, indicating that the additional vertical shift does not 
play a significant role in this case. In fact, the increased oscillations in the peel stress distribution 
caused by the vertical shift may not be desirable in some situations. Alternatively, when the 
fibers are vertically shifted without decreased horizontal spacing, the peel stress is reduced by a 
modest 9%. Thus it appears that the major contribution to the reduction of the free-edge peel 
stress comes from functional grading in the horizontal direction. 

It is interesting to relate the reduction in the peel stress produced by the functional grading 
of fiber architecture to the reduction in the transverse thermal expansion mismatch between 
B/Ep and Ti plies with increasing fiber volume fraction illustrated in Figure 4. As observed in 
the figure, the transverse thermal expansion mismatch decreases rapidly with increasing fiber 
volume content for fiber volume fractions greater than approximately 0.05. As observed in Table 
4, the local fiber volume fractions in the first three vertical columns of cells in the uniformly 
spaced configurations are 0.55, 0.50 and 0.50, compared to 0.63, 0.60 and 0.56 in the linearly 
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spaced configurations. Figure 4 indicates that the transverse thermal expansion coefficients in a 
B/Ep ply with fiber volume fractions of 0.50, 0.55 and 0.63 are 32.50, 29.52 and 24.91 
(10“V > C), respectively. The relatively modest decrease in the local thermal expansion mismatch 
(approximately 14%) in the immediate vicinity of the free edge (i.e., the first three fibers) 
achieved by a correspondingly modest increase in the local fiber volume fraction (approximately 
14%) thus produces a substantial reduction in the peel stress at the free edge (approximately 
25%). This obviously has technologically significant implications. 

4.0 SUMMARY AND CONCLUSIONS 

A previously developed theory for the elastic response of metal matrix composites with a 
finite number of uniformly or nonuniformly spaced inclusions or fibers in the thickness direction 
subjected to a through-thickness thermal gradient has been extended herein to enable analysis of 
mate rial architectures characterized by nonuniform fiber spacing in two directions. In this new 
approach, the microstructural and macrostructural details are explicitly coupled when solving 
the thermomechanical boundary-value problem. Coupling of the local and global analyses 
allows one to rationally analyze the response of polymeric and metal matrix composites such as 
B/Ep, B/Al and SiC/TiAl that contain relatively few through-thickness fibers, as well as so- 
called functionally graded materials with continuously changing properties due to nonuniform 
fiber spacing or the presence of several phases. For this class of emerging composites, it is diffi- 
cult, if not impossible, to define the representative volume element (RVE) used in the traditional 
micromechanical analyses of macroscopically homogeneous composites (Hill, 1963). 

The extension of the theory to material architectures functionally graded in two directions 
makes possible the analysis of laminated composite plates with finite dimensions in one plane 
subjected to combined thermo-mechanical loading. In particular, the technologically important 
interlaminar stress fields in laminated composite plates in the vicinity of the free edge can be 
analyzed with the outlined approach. The presented comparison of interlaminar stresses in a 
symmetrically laminated B/Ep-Ti plate subjected to a uniform temperature change generated 
using the new coupled approach and the finite-element method based on homogenized properties 
of the individual plies demonstrates the capability of the proposed theory to capture the large 
stress gradients at the free edge. Explicit incorporation of the microstructure of the B/Ep plies in 
the interlaminar stress calculations using the coupled theory produces maximum value for the 
peel stress at the free edge that is approximately 35% higher than the corresponding value based 
on the homogenized B/Ep ply properties (Figure 8). This requires explicit consideration of at 
least seven columns of fibers localized at the free edge. Microstructures with fewer fibers pro- 
duce stress fields that do not adequately reflect the actual stresses. Thus calculation of the 
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technologically important interlaminar peel stress at the free dege of a laminate based on the 
homogenized continuum approach substantially underestimates the actual stress fields, leading 
to unsafe designs. 

In order to reduce the high interlaminar peel stress at the free edge, the local thermal expan- 
sion mismatch between the B/Ep and Ti plies must be reduced. This can be achieved by altering 
the microstructure of the B/Ep plies at the free edge through functional grading of the boron 
fiber distribution. By decreasing the fiber spacing with decreasing distance from the free edge 
the local fiber volume fraction in the B/Ep plies is increased, thus decreasing the transverse ther- 
mal expansion mismatch as observed in Figure 4. Substantial reductions in the free-edge peel 
stress can be achieved by a relatively modest increase in the local fiber volume fraction directly 
at the free edge (see Figure 13 and Table 4). 

It should be emphasized that the new theory presented herein makes possible the investiga- 
tion of the effects of fiber distribution and fiber shape in functionally graded composites. Recent 
investigation of these effects in doubly-periodic composites was conducted by Arnold et al. 
(1994) using the concept of a repeating unit cell. Such effects can now be investigated in the 
presence of a continuously changing microstructure in conjunction with the coupling of the 
micro and macro-structural response. 

Finally, in previous investigations the authors have demonstrated that inclusion of inelastic 
effects in the analysis of functionally graded materials may be important in the presence of large 
temperature gradients in applications involving one-dimensional version of the higher-order 
theory (cf. Aboudi et al., 1994c, Pindera et al., 1995). Others have demonstrated the importance 
of inelastic effects within the context of the free-edge problem. For example, Williamson et al. 
(1993) and Drake et al. (1993) employed the finite-element method to study the residual stresses 
that develop at graded ceramic-metal interfaces joining cylindrical bodies made of metallic and 
ceramic materials. The gradation was modeled using a series of perfectly bonded cylindrical 
layers, with each layer having slightly different properties. Their results demonstrate the impor- 
tance of plasticity effects in the analysis of graded and non-graded interfaces. The authors also 
showed that in some cases optimization of the microstructure of graded layers is required to 
achieve reduction in certain critical stress components that control interfacial failure. Along 
similar lines, Suresh et al. (1994) studied the response of elastoplastic bi-material strips sub- 
jected to cyclic temperature variations. Closed-form solutions were derived using simple beam 
and plate theories to analyze the stresses and curvature that develop in a bi-material configura- 
tion for the considered thermal loading. In addition, finite-element formulation was employed to 
capture other features not included in the simple analytical models. In particular, the authors 
showed that the plastic flow along the interface separating the two materials at the free edge can 
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be modified substantially by altering the constraints at the edges of the strip. Therefore, in subse- 
quent investigations inelasticity effects will be incorporated into the theoretical framework of 

HOTFGM-2D. 
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7.0 APPENDICES 


7.1 Appendix A: Thermal Analysis 


A.1 Heat Conduction Equations 

Multiplying equation (1) by (it) / (x£ ) m (*3 ) n > where /, tn, n = 0, 1, or 2 with l +m + n ^ 2, and 
integrating by parts, using the temperature expansion given in equation (9), the following equations are obtained: 


+ L $$0) + I»,0) =0 < A1 > 

^.€)-eStt«=0 (A2) 

M#i)-(2^,0)=0 (A3) 

y4 [ 3L^, 0 ) + 0) + *M,0) 1 - 2fiWU = 0 (A4) 

y/i£? )2 [ M# o) + 3L^, 0) + i».o) 1 - 20$fto) = 0 (A5) 

y/<T )2 [ M# 0) + ^,0) + 3L^Uo> ] - 20^,1) = 0 (A6) 

where Q.\^n >n ) has been defined previously in equation (10), and 
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(A8) 


(A9) 


d hW } r 

In the above: n = 0 or 1; 9^ } (±y-). 9 ^(±-y) denote the interfacial fluxes at 

3q a) = +— d a , xf ] = ±— h$\ xf = ±—l ( y\ respectively; and v^ r) = ^a/i^ ^ r) is the volume of the subcell 
2 2 2 

(afty) in the (p,q,r)th cell. 

Equations (Al) through (A6) provide relations between the zeroth-order and first order heat fluxes 0f“$£ n ) 
and the interfacial fluxes Explicit expressions for the interfacial fluxes L\f^ n ) given solely in terms of 



Gi$/n n) 316 obtained through the following sequence of manipulations, noting that equations (A2) and (A3) already 
provide direct relations between and £$$,0)- and and First, substituting equation (Al) 

into (A5) and (A6), respectively, gives the following direct expressions for and L$)$ 0) . 

M#0)=12e^,0)^{f )2 (A10) 

Utti.O) = 120$$.i)4 )2 (All) 

Then upon substitution of equations (A10) and (All) into (Al) we obtain the following expression for 

Lffiti, 0) = -12(e^,0)^^ )2 + <2$$,lA r)2 ) (A 12) 

Equations (A 10) through (A 12) will be used to reduce the heat conduction and heat flux continuity equations to 
expressions involving only the heat flux quantities These can subsequently be expressed in terms of the 

fundamental unknown coefficients 7*$$, appearing in the temperature expansion given by equations (9), using 
equations (11) through (15). 

A.2 Heat Flux Continuity Equations 

The heat flux continuity conditions (3)-(4) are imposed on an average basis at each subcell and cell interface. 
The heat flux condition in the x \ direction is obtained usingequation (3a) in (A7), yielding 


(AB) 

We note that equation (4a) is identically satisfied for a material that is periodic in the x\ direction by the chosen 
temperature field representation. 

Prior to imposing the continuity conditions in the FG x 2 and x 3 directions, let us define intermediate quanti- 
ties and as follows: 


ie* 1 *>*-*?'» i 

«**• i lp '"’ r) = ? r> & + i 


(p,q,r) 

Xi =-hf ) /2 
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jtl) 




(A14) 
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These q uan tities will simplify the algebra associated with application of the heat flux continuity requirement on an 
average basis in the x 2 FG direction. Then substituting equations (3b) and (4b) into the above definitions we have, 
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(A16) 
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Adding and subtracting equal quantities to and from equations (A16) and (A17) it can easily be verified that 

2/$“ l7) I iP ' qr) = + g^ t2 ^] (p ' q ' r ' > - {ff 2 ^ +g^ ot 2 7 )](P.?-l.r) (A 18) 

2g% xly) I (P ' qr) = [-fti 2 ^ +g < *2* b(i ] (p,9,r) + + g£ x2Y) ] <p ' 9 ~ , ’ r) (Ai9) 

Then using equations (A18) and (A19) in (A 8 ), we obtain the following heat flux continuity conditions for the x 2 
FG direction: 

i -h ,m°> * m.o> - ^ - 1 L ®-»> + t z -® 8 -°> 1<p ' , ' u> - 0 <A20> 

[ -iffl. 0 ) + yiffl.o, - %-iM 0) ^- r) + jl + = 0 <“» 

Similarly, the heat flux continuity conditions (3c) and (4c) in the remaining x 3 FG direction provide 


and 


+um» - 4 «» =° <A22) 



Equations (A10) through (A12), together with equations (A2) and (A3), will be employed in reducing the heat 
flux continuity conditions (A13), (A20) through (A23) to expressions involving only the volume-averaged zeroth- 
order and first-order heat flux quantities 


A.3 Reduction of Heat Conduction and Heat Flux Continuity Equations 

Substituting equation (Al) into (A4), and using equations (A10) and (All), reduces the volume-averaged heat 
conduction equations to a set of 8 equations given by equation (16). Next, using the expression for L\ (§Tfj,o) given 
by equation (A12), and the expression for L^?,o) given b y equation (A2), and given by equation (A3), in 

the continuity relations (A13), (A20) through (A23), we obtain equations (17) and (21). 
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A.4 Thermal Continuity Equations 


As in the case of the heat flux field, the thermal continuity conditions (5) - (6) are imposed on an average 
basis at each subcell and cell interface. Thus substituting equation (9) into equations (5a) - (5c), respectively, we 
obtain at each subcell interface the conditions (22), (23) and (25). Furthermore, substituting equation (9) into (6b) - 
(6c) to ensure continuity of temperature between neighboring cells in the FG directions, we obtain the conditions 
(24) and (26). 

7.2 Appendix B: Mechanical Analysis 
B . 1 Equations of Equilibrium 

Let us multiply equation (28) by (xf *) 1 (x^)™ (S^V. where again /, m, n = 0, 1, or 2 with 
/ + m + n < 2. Integrating the resulting equations by parts, and using the displacement expansions (37), we obtain 
the equations of equilibrium in the subcell region (aJ3y) in the form: 


4;^?o,0) + 45to,0) + ^tj$3?o,o) = 0 

(Bl) 
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(B2) 
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~l2^a [ 3/\“^o,0) + J %j$!o,o) + ^^3?o.o) ] ~ 25^^o,o) = 0 

(B5) 

[ A/f&?o,0) + 3 ^“^o,0) + ^Is^o.o) ] ~ 25^3?i,o) = 0 

(B6) 

^•4 r)2 [A® 0.0) + 4$?0.0) +4^0,0)1-25^0,1) =0 
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where j = 2,3 in equations (Bl) and (B3) through (B7), while in equation (B2) j = 1. In the above equations, 
S [ has been defined previously in equation (38), and 
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d a /2 hf/1 

= J } [oSf«(V)+(-i)' t '^r(44' ) )]^. a w (BIO) 
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In the above: n = 0 or 1; G\f^(±^d a ), 0^(±j4 r) ), stand for the interfacial stresses at 

3c ( | a) =+^ a ,lf =±-^h$\xf =±y4 r) , respectively. 

Equations (Bl) through (B7) provide relations between the zeroth-order and first order, volume-averaged 
stresses and the interfacial tractions /\“^ 0 ,0). 4j$?n,0) and ,«)• Direct "one-to-one" relations are 

obtained through the following sequence of manipulations, noting that equations (B2) through (B4) already provide 
direct relations between Snfj^o) and ^Vipi^o.O) • and ar *d 1 ), j = 2,3. First, 

substituting equation (Bl) into equations (B6) and (B7), respectively, gives direct expressions for -Ayp^o.O) and 

*»,0). 

4^0,0) = 12S^> 1 ,0)^ )2 (BID 

M^o, 0 ) = 125^o.t)4 )2 < b12 ) 

Then, upon substitution of equations (Bl 1) and (B12) into (Bl), we have the following expression for /$$o,0) : 

/$3?0,0) =-12 (S», 0) ^ )2 +^o,i)4 )2 ) (BID 

Equations (B 11) through (B13) will be used to reduce the equilibrium equations and traction continuity equations to 
expressions involving only the zeroth-order and first-order, volume-averaged stress quantities These can 

subsequently be expressed in terms of the fundamental unknown coefficients appearing in the displacement 

field expansion given by equations (37), using equations (39) through (46). 

B.2 Traction Continuity Conditions 

The traction continuity conditions, equations (31) - (32), are imposed on an average basis at the subcell and 
cell interfaces. These conditions imply existence of certain relationships between the surface integrals of the interfa- 
cial traction components defined by equations (B8) - (BIO). The normal traction continuity condition in the X\ 
periodic direction is obtained using equation (31a) in (B8) with j = 1. This yields 

< B14 > 
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We note that equation (32a) is identically satisfied for j = 1. 

For the shear traction continuity conditions in the Xj direction we obtain from equation (31a) and (B8) with j 

= 2,3 


U,/Wo.® +d 2 i\fS M (BIS) 

Again, equation (32a) is identically satisfied for j = 2,3. 

To assist in establishing the continuity relations in the FG directions, let us define two new quantities 
and as follows 
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Substituting equations (31b) and (32b) into the above definitions, we obtain, respectively: 
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By addition and subtraction of equal quantities to and from equations (B18) and (B19) it can easily be verified that 


2 F^“ lr) I {p ' q ' r) - [_J irfe^T) + (J^2y)](P.9.r) _ + (jl« 2 Td](P.9-l.r) (B20) 

2G^“ lY) | (P <? r) _ _ /r^2Y)](P.9.r) + [G^ + F%f r ' ) ] (p ’ q ~ 1 ' r) (®21) 

Then employing equations (B20) and (B21) in equation (B9) with j = 2,3 we obtain the corresponding relations: 


[ -h i/gWo.® + - jkiAffio.O) - 1 JffltM + y*2^Jo,o,F' , ‘ l ' r) = 0 (B22) 
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Similarly, using equations (31c) and (32c) with j = 2,3 we have 
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As a result of the above manipulations, 44 relations, given by equations (B14) - (B15) and (B22) - (B25), arise 
from the traction continuity conditions between subcells and between neighboring cells. These equations, in con- 
junction with equations (B2) through (B4), and equations (B 11) through (B13), will be employed in reducing the 
equilibrium and traction continuity equations to expressions involving only volume-averaged zeroth-order and first- 
order stresses 


B.3 Reduction of Equilibrium and Traction Continuity Equations 

Substituting equation (B 1) into equation (B5) and using equations (Bll) and (B12), reduces the volume- 
averaged equilibrium equations to a set of 16 equations given by equation (47). Next, combining the expressions for 
itfflfoo) and n*0,0). provided by equations (B2) and (B13), and the continuity relations (B14) - (B15), respec- 
tively, we obtain the twelve equations given by (48) to (50). Continuing, if we substitute equation (B3) and (Bll) 
into equations (B22) - (B23) directly, we obtain equations ((51) and (52), respectively. Finally, combining equations 
(B4) and (B12), and equations (B24) and (B25), yields equations (53) and (54), respectively. 


B.4 Displacement Continuity Conditions 

The displacement continuity conditions, i.e., equations (33) - (34), are now imposed on an average basis at the 
interfaces. This is accomplished by first substituting equation (37) into equation (33a), yielding equations (55) - 
(57), then into equations (33b) - (34b), yielding equations (58) - (61), followed by equations (33c) - (34c), yielding 
equations (62) and (65). Consequently, equations (55) - (65) provide 44 relations which must be imposed to guaran- 
tee the continuity of the displacements between the subcells and between neighboring cells. 
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Table 1. Material properties of the boron/epoxy plies (v/ = 0.50) and titanium sheets. 


Material 

E* a (GPa) E* t (GPa) 

* 

Va 

a* (10 _6 m/m/°C) 

oct-OO - 6 m/m/°C) 

B/Ep 

206.8 18.6 

0.21 

4.5 

30.6 

Ti 

118.6 118.6 

0.34 

8.5 

8.5 


Note: Subindices A and T denote axial and transverse quantities, respectively. 


Table 2. Material properties of boron fiber and epoxy matrix. 


Material 

E (GPa) 

V 

a(10 _6 m/m/°C) 

Boron fiber 

417.0 

0.13 

6.0 

Epoxy matrix 

5.24 

0.35 

50.0 


Note: E and v denote the Young’s modulus and Poisson’s ratio, respectively, and a is the coefficient of 
thermal expansion. 


Table 3. Center-to-center distances between fibers in the cells at the free edge in the linearly spaced confi- 
gurations. 


Adjacent fibers from the free edge 

Center-to-center distance between adjacent fibers (|i.m) 

Free edge - 1 

97.36 

1-2 

201.08 

2-3 

215.90 

3-4 

230.70 

4-5 

245.50 

5-6 

260.35 

6-7 

275.20 

7-8 

290.00 

8 - homogenized B/Ep 

215.90 
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Table 4. Local fiber volume fractions in the cells at the free edge. 


Cell number from free edge 

v f (uniform fiber spacing) 

v f (linearly increasing fiber spacing) 

1 

0.55 

0.63 

2 

0.50 

0.60 

3 

0.50 

0.56 

4 

0.50 

0.52 

5 

0.50 

0.49 

6 

0.50 

0.46 

7 

0.50 

0.44 

8 

0.43 

0.34 


Note: The fiber volume fraction in the 1st and 8th cells of the uniformly spaced configuration differs from 
0.50 due to horizontal shift of the fiber architecture towards the free edge as explained in the text. 
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Figure 2 - (a) Three-dimensional schematic of a composite functionally graded in the x 2 
and x 3 directions showing the dimensions of the basic building block (b) of the composite. 
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Figure 4 - Effective normal and transverse thermal expansion coefficients of B/Ep as a 
function of the fiber volume fraction. Also shown in the figure is the CTE of the isotropic 
titanium sheet. 
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Figure 5 - Normal and shear stress distributions in the titanium ply at the interface 
separating B/Ep and Ti plies due to a temperature change of AT* = - 154.45°C. Comparison 
between the predictions of HOTFGM-2D and FE analysis. 
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Figure 5 (cont’d) - Normal and shear stress distributions in the titanium ply at the interface 
separating B/Ep and Ti plies due to a temperature change of A T=- 154.45°C. Comparison 
between the predictions of HOTFGM-2D and FE analysis. 
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Figure 6 - Schematic of eight configurations of the [ B/Ep-Ti ] s laminate showing increas 
ingly refined microstructure at the free edge. 
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Figure 7 - Normal stress 022 distributions in the titanium ply at the interface separating 
B/Ep and Ti plies due to a temperature change of A7" = - 154.45°C for the eight configura- 
tions of Figure 6. 
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Figure 7 (cont’d) - Normal stress 022 distributions in the titanium ply at the interface 
separating B/Ep and Ti plies due to a temperature change of AT = - 154.45°C for the eight 
configurations of Figure 6. 








Figure 8 - Normalized maximum stress 022 in the tita ni u m ply at the free-edge of the lam- 
inate at the interface separating B/Ep and Ti plies due to a temperature change of AT — 
154.45°C as a function of the refined microstructure for the eight configurations of Figure 
6. The normalization is performed with respect to the corresponding maximum stress 
obtained from the analysis based on homogenized B/Ep properties. 
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Figure 9 - Effect of removing fibers near the free -edge of the laminate on the normal stress 
022 distribution in the titanium ply at the interface separating B/Ep and Ti plies due to a 
temperature change of AT = - 154.45°C. 
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fiber #2 removed 


■ fibers #2 and #3 removed 


Figure 10 - Ma ximum stress 022 m the titanium ply at the free-edge of the laminate at the 
interface separating B/Ep and Ti plies due to a temperature change of AT = - 154.45°C 
obtained from distributions of Figure 9, normalized with respect to the corresponding max- 
imum stress obtained from the configuration with uniformly-spaced fibers. 
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Figure 11 - Schematic of different fiber distributions in the B/Ep ply near the free edge 
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^ uniform spacing - vertical shift 
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Figure 13 - Maximum stress 022 in the titanium ply at the free-edge of the laminate at the 
interface separating B/Ep and Ti plies due to a temperature change of AT = - 154.45°C 
ob tain ed from distributions of Figure 12, normalized with respect to the corresponding 
maximum stress obtained from the configuration with uniformly-spaced fibers. 
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